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The author analyzes quantum computation with the hybrid qubit (HQ) that is encoded using 
the three-electron configuration of a double quantum dot. All gate operations are controlled with 
electric signals, while the qubit remains at an optimal operation point that is insensitive to noise. An 
effective single-qubit description is derived, and two-qubit interactions are suggested using Coulomb 
and exchange interactions. Universal quantum control is described and numerically simulated using 
realistic parameters for HQs in Si and GaAs. High-fidelity quantum computing at the threshold of 
quantum error correction is possible if the Coulomb interactions between the HQs stay weak. 


I. INTRODUCTION 

Spin qubit quantum computers are pro mising plat¬ 
forms to achieve quantum computation.^^^ The electron 
spin of a gate-defined quantum dot (QD) naturally de¬ 
fines a two-level quantu m system tha t en codes one bit 
of quantum information.l^I^ Also a duo,l^ a trio,^ and a 
quartet!^ of such singly occupied QDs have been pro¬ 
posed as realizations of quantum bits. The so-called hy¬ 
brid qubit (HQ) is encoded using a double QD (DQD) 
that is occupied with three electrons.!^ The name HQ 
was introduced because the HQ is either a spin or a 
charge qubit depending on its operation principle. A spin 
qubit is well protected from charge noise, but typical spin 
qubits cannot be operated all-electrically. A charge qubit 
can be manipulated with electric fields, but it is therefore 
also susceptible to charge noise. 

The HQ can be operated exclusively with electric sig¬ 
nals similar to the triple QD (TOD) qubit that is coded 
using three singly occupied QDs.EIS Because subnanosec¬ 
ond con trols of the electrostatic potentials of QDs can be 
realized,electric manipulations of QD qubits are al¬ 
ways more favorable than magnetic manipulations. Even 
though sin gle sp ins can also be controlled with magnetic 
field pulsesjg difficult to selectively control a single 
spin with magnetic fields when this spin is in the vicin¬ 
ity of other spins .1^^ Electric fields can also modify single 
spins indirectly (e.g. via spin-orbit interactions or mag¬ 
netic field gradients), but high-fidelity gates still remain 
challenging .EHIZl 

Most conveniently, the HQ is manipulated by the 
transfer between the (u-qDiiJ^qDj) = (1)2) and (2,1) 
configurations similar to a charge qubit 

and are the electron numbers at QDi and QD 2 . 

These manipulations require precise, subnanosecond 
pulses of electric signals with large amplitudes. Consid¬ 
erable advances for the HQ were possible using an asym¬ 
metric DQD configuration .E 2 HIII Por a doubly occupied 
QD, the addition of the second electron usually requires 
higher energy in the triple configuration compared to the 
singlet configuration. DQDs were constructed, where this 
singlet-triplet energy difference is small for one of the 


QDs, but it is much larger for the other one. 

This paper studies DQDs with asymmetric addition 
energies in the layout of Fig. QD 2 always has a small 
singlet-triplet energy difference, but the singlet-triplet 
energy difference difference is large for QDi. It will be 
shown that the HQ has a small energy difference in the 
(1,2) configuration. When approaching the (2,1) config¬ 
uration, the qubit states pass through two avoided level 
crossings. It h as bee n described theoreticallj® and shown 
experimentallj®^^ that these two anticrossings are suf¬ 
ficient for the single-qubit control of the HQ with the 
appropriate tuning protocols through the anticrossings. 

It is also possible to operate the HQ in the vicinity 
of an anticrossing, where all the qubit operations are 
realized with microwave pulses of small amplitudes.!^ 
An anticrossing is a sweet spot in the energy diagram 
where low-frequency noise in the control parameter does 
not dephase the qubit. Such sweet-spot operations 
have improved the coherence times for superconducting 
qubits.f^^t^ZlNote that microwave gates have also been re¬ 
alized for ot her spin qubit encodings, as for DQDd^ and 
TQDs .1221331 

For HQs, the entangling operations remain more chal¬ 
lenging than the single-qubit gates. Ref. jH] proposed 
that strong electrostatic couplings can be used to entan¬ 
gle two HQs. Also the rapid transfer of electrons between 
the HQs has been studied.l^ In this paper, I show that 
the HQ can always be operated at its sweet spot while 
weak couplings between HQs enable two-qubit gates. I 
quantify the inter-qubit couplings with Coulomb and ex¬ 
change interactions, and show that both interactions en¬ 
able two-qubit gates. 

The main findings of this paper are explicit manipu¬ 
lation protocols for HQs of high fidelities, while all the 
HQs are operated at their sweet spots. I simulate the 
gate operations with typical qubit parameters and show 
that these gates tolerate the dominant noise sources of 
HQs. At the sweet spot, the HQ keeps enough character 
of a spin qubit to be protected from charge noise, but it is 
already sufficiently close to a charge qubit to be rapidly 
tunable. To employ these gates experimentally, one has 
to build setups that have weak Coulomb interactions be- 
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Figure 1. Setup of two HQs. Each HQ is coded using the 
three electron configuration of a DQD. The superscripts L 
and R label the positions of the HQs; the subscripts 1 and 2 
label the QDs. 


tween the HQs at their sweet spots. 

The organization of the paper is as follows. Sec. |TT] in¬ 
troduces the description of the HQ, and it specifies the 
magnitudes of the relevant parameters. The single-qubit 
control, the readout, and the initialization are also de¬ 
scribed. Sec. lHII derives the effective Hamiltonians for the 
couplings between HQs that are generated from Coulomb 
and exchange interactions. Sec. |IV| specifies two ap¬ 
proaches to realize universal control of the HQ. Not only 
entangling operations are discussed, but it is also shown 
that single-qubit gates are possible. Sec. |V] discusses the 
influence of electric and magnetic noise for HQs that are 
operated at their sweet spots. Sec. |VI| summarizes the 
findings. 


II. SINGLE-QUBIT DESCRIPTION 


The HQ is coded using a DQD that is occupied with 
three electrons, as shown in Fig. The qubit manipula¬ 
tions require a much larger singlet-triplet energy splitting 
of a doubly occupied QD for QDi compared to QD 2 . 
I consider only the subspace enc oding in the S = 1/2, 
Sz = 1/2 spin configuration.lSl^ Even though a global 
magnetic field is not required to realize single-qubit gates, 
it will be necessary to apply magnetic fields to initial¬ 
ize a HQ to S' = 1/2, = 1/2.1^ QDi can be used 

for the initialization of a singlet because a two-electron 
configuration quickly equilibrates to a singlet under ther¬ 
mal relaxation. In the presence of an external magnetic 
field, also the spin-up configuration at QD 2 can be ini¬ 
tialized because it has lower energy than the spin-down 
configuration .1^ 

The HQ is described in the basis of the lowest states 
with S = 1/2, Sz = 1/2. One considers |x) = |S f) in 
the (2,1) configuration, and the states 11) cx \ T_|_) — 

ItTo) and | 0) = jfS) in the (1,2) configuration. The 
singlet |S) oc Ifi) — lit) and the triplets |T+) = |tt)) 
I To) (X Itt) + lit)) and |T_) = |tt) describe the two- 
electron configurations at a QD. The three-electron states 
with S = 1/2, Sz = 1/2 are defined using the standard 
spin addition rules .1^ The Hamiltonian for the HQ in the 


(2,i)<— -*aa) 



Figure 2. Energy diagram of the relevant states to describe 
the HQ in the S = 1/2, Sz = 1/2 spin configuration of a three- 
electron DQD, according to Eq. 0. for 2Ai/h = 14.5 GHz, 
2Ao/h = 5.2 GHz, and 5/h = 12.1 GHz. e models the energy 
detuning between the (2,1) and (1,2) charge configurations. 
The HQ is coded using the two states with the lowest energies. 
The HQ is in (1,2) for e —>■ cxa with the excited state | 1) oc 
y/2 1 4 , T+) — I h To) and the ground state | 0) = 11 S). Only 
the ground state | x) = | S t) is in (2,1) for e —^ —00, while the 
excited state remains in (1,2). There are two state couplings 
(I 1) -f-)- I x) and I 0) | x)) at e ~ 0. The gray line shows the 

S = 3/2, Sz = 1/2 state (called | y)) in (1, 2) that is uncoupled 
to (2,1). All gate operations for the HQ are realized at the 
sweet spot e*, where the qubit is insensitive to small variations 
in e. The dashed lines approximate the qubit levels in the 
vicinity of e, [see Eq. §]• The parameters Ai and Ao are 
defined in the text. 


basis {|x) , 11) , 10) } is 


/ e -Ai Ao\ 

'H{|x).|i).|o)} = -Ai 8 0 . (1) 

V Ao 0 0 y 

The Coulomb repulsion is high for a doubly occupied 
QD in the triplet configuration, which raises the energy of 
11) compared to |0) by 5. |x) and 11) (|x) and |0)) are 
coupled via the real parameter Ai > 0 (Aq > 0), which 
is called the tunnel coupling in the following, e is the 
detuning parameter between the (2,1) and (1, 2) config¬ 
urations. The electron configuration in (2,1) [(1,2)] has 
lower energy for e < 0 [e > 0], and (2,1) and (1,2) have 
equal energies at e = 0. 

Eq. 0 describes a three level system with two anti¬ 
crossings. Fig.|^shows the energy levels for typical qubit 
parameters in the transition region between (2,1) and 
(1,2). The two states with the lowest energies encode 
the HQ. For e <C 0, only the ground state is in (2,1), 
while the first excited state is in (1,2). This configura¬ 
tion is ideal for the readout of the HQ because a nearby 
charge detector can distinguish the two states of the HQ, 
similar to the readout of a singlet-triplet qubit por 
e 0, the HQ is in (1,2) and the energy difference be¬ 
tween the qubit states is unchanged when e is slightly 
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modi fied. This regime is highly protected from charge 
noise .1221401 Note that there is additionally one S = 3/2, 
S 2 = 1/2 state in the ( 1 , 2 ) configuration (called | y)) that 
is uncoupled to the ( 2 , 1 ) configuration. 

I consider a HQ with similar magnitudes of the singlet- 
triplet energy splitting at QD 2 S and the tunnel couplings 
Ai and Ag: S > AijAg. The singlet-triplet energy 
splitting at QDi is much larger than S such that dou¬ 
bly occupied triplets at QDi are neglected. Note that 
such HQs have been realized using Si DQDs,l2^^^ where 
the singlet-triplet energy difference depends on the or¬ 
bital and the valley energy splittings. The orbital energy 
splitting is determined by the confining potential of a 
QD, while the valley energy splitting is determined by 
the potential landscape at the atomic scale. A recent ex¬ 
periment showed that the valley splitting of a Si QD can 
be controlled electrically.SH Note that HQs can also be 
build using GaAs QDs, but these setups necessarily re¬ 
quire asymmetric DQD potentials to operate them with 
the equivalent tuning pulses. In this case, the smaller QD 
has a large singlet-triplet splitting; the larger or elliptic 
QD has a small singlet-triplet splitting.ESHlll The param¬ 
eters of previous experiments with Si HQs are used in the 
following: 2Ai//i = 14.5 GHz, 2Ao/h = 5.2 GHz, and 
d/h = 12.1 GHz (cf. Refs. |19l2nl21| l. 

Ai in Eq. 0 is caused by the coupling to an ex¬ 
ited orbital, which is larger than the ground state or¬ 
bital that determines Ag. Tunnel couplings are expo¬ 
nentially suppressed with the distance between localized 
orbitals. Therefore, Ai is significantly larger than Ag, 
and it is meaningful to analyze Eq. Q in the eigen- 
basis of the larger anticrossing. The states are ro- 


1 = 


tated to |x) = cos ( 0 e/ 2 ) |x) -f sin {0^/2) \ 1 ) and 
— sin ( 0 e/ 2 ) |x) -l-cos {0^/2) \ 1 ), with sin [9^) = — 2 Ai/C/e, 


cos {9,) = (e - A)/C/„ and U, = sj -f {e-Sf. 
This basis rotation modifies Eq. 0 to 
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The coupling of | x) to | 1 ^ , 10) | is neglected in 

Eq. ([^ to describe the low-energy subspace that defines 
the HQ. It will be shown that the latter subspace has 
an anticrossing at e* = Af/S, while U^, ^ Ag. The 

two-level system | 1 ^ , | 0 ) | is approximated by: 


with Ag = 2 Aosin( 6 »*/ 2 ), A{e) = 2 (^ 2 ^^) ~ 

and 0, = . To derive Eq. 0 , the dependency 

of sin {9f^/2) on e is neglected because this term varies 
slowly in the vicinity of e*. Then the qubit’s energy 


{|T),|0)} 


splitting 4 ) + [^0 sin ( 0 *)] within 

is minimized with respect to e, giving the first term 
in Eq. 0, with Ag = 2Agsin(0»/2) for e» « 

A\/5. The second term in Eq. ^ is obtained from 


de [Ag sin {9^)Y 

The basis in Eq. ^ has been rotated to 




T> 
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1 
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|x) +COS 
|x) +COS 


|1) +|0) 

|1) -|0) 


(4) 

(5) 


with the definitions of the Pauli operators = |l) (l| — 
|0) (0| and ax = |l) (0| -I- |0) (l|. Note that 11) and 
I 0 ) have finite contributions in ( 2 , 1 ) and ( 1 , 2 ). 0 * « 
1 . 347 r is the mixing angl^ with cos^ (^*/2) ~ 0.26 and 
sin^ (0*/2) « 0.74. Fig. proofs that the effective two 
level system in Eq. ([^ describes the HQ close to e». The 
leakage state |x) is raised hy Ai = 5 + A\/5. Ag = 
2Agsin(0*/2) is the energy difference between 11) and 
| 0 ). 

Eq. ^ permits the usual Rabi control of the qubit 
with microwave drives of the detuning parameter gHSHSl 
If e is driven with small amplitudes around e* [A (e) —> 
Acos (27rDt//i -I- (/))], then all possible single-qubit oper¬ 
ations can be realized for D = Ag when the phase is 
varied. Eq. ([^ gives in the rotating frame with the 

static Hamiltonian 


y^rwa. _ ^ _ ^ 0 ^ 

Eq. (|^ uses implicitly a rotating wave approximation, 
which is valid for Ag ^ A. 


III. TWO-QUBIT INTERACTIONS 

This section describes the interactions between two 
HQs. The superscripts L and R label the positions of 
the HQs. QD 2 and QD? are the neighboring QDs from 
HQL and HQR (cf. Fig.0. 


A. Capacitive Coupling 

The charge configurations of two HQs couple capaci- 
tively via the Coulomb interaction. The dominant con¬ 
tribution is determined by the electron configurations at 
QD? (riQigL) and QD? (figof) according to 

C = K ^QD^^QDf ■ C^) 

The magnitude of the coupling parameter k depends on 
the layout of the experiment and the QD material, n can 
be large; e.g., one can approximate its magnitude by the 
Coulomb interaction of two electric point charges that 
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are 250 nm apart, giving k 


500 /xeV for 


47reo^r 250 nm 

the dielectric constant Cr = 11.7 of silicon .1^ This value 
agrees with the approximation of a few tenth of meV in 
Ref. [48] . 

The entangling operations that follow require much 
smaller k than the naive estimate given above. The envi¬ 
ronment around the and the metallic gates partly 

screen the capacitive couplings between the HQs, and 
thus they reduce k. Also the layout of the four QDs can 
be be designed such that the Coulomb interactions be¬ 
tween neighboring HQs are lowered.!^ Finally, bringing 
QD 2 and further apart lowers k. 

For two HQs operated at their sweet spots, 
and e?, the projection of Eq. (|7| to the subspace 

the two- 


|HQ) 


IHQ^HQ**) 




qubit interaction 


oLqRX 


gives 


C 2 Q — 


( 8 ) 


with Ac = -f sin^ (0L/2)sin^ (6»^/2). Eq. ([§ lowers 


TLtR 


0 "T^) + 
0 "T^\ - 


^Lt^R 


tLtR 




and 


j and 


the energies of the configurations 
compared to 

1""0“^ . The first states have large weights in 

( 2 , 1 , 2 , 1 ) or ( 1 , 2 , 1 , 2 ), while the latter states have large 
weights in ( 1 , 2 , 2 , 1 ) or ( 2 , 1 , 1 , 2 ). 

Additionally to Eq. (|^, there are also the single-qubit 
interactions 


CiQ =Xc<j]^ + 


rR R 

'X 7 


(9) 


with = -| sin^ ( 6 »l/ 2 ) sin^ ( 6 '?/ 2 ) and = |[1 -h 
cos^ {9l/2)]s\T? {ef/2). 


B. Exchange Coupling 

Encoded qubits in the multielectron regime have a 
large number of spin states that are not part of the com¬ 
putational subspace.l^ As a consequence, the time evolu¬ 
tion out of the computatio nal su bspace, which is called 
leakage, must be considered .1^*^ Similar to TQDs, inter¬ 
qubit tunnel couplings between HQs can cause leakage 
because the subspac es o f different local spin quantum 
numbers are coupled.It will be shown that the sys¬ 
tem of two HQs at their sweet spots have large energy 
separations between leakage and computational states, 
which suppresses leakage events (cf. similar approaches 
for TQDs in Refs. |53l54| h 

In the six-electron configuration, only states with total 
Sz = 1 are considered because the spin directions are con¬ 
served during tunneling events. Tab. [^ summarizes the 
spin states with the lowest energies. The computational 
subspace | 1 ^ 1 ^^ , , 0 ^ 0 ^^ | and the 

leakage states (called | a), \j3), 1 7 ), 115), and | v )) are con¬ 
sidered. Only the states with 5” = 1/2 are able to couple 





Figure 3. States with the lowest energies for a HQ and 
for the combined system of two identical HQs in the ab¬ 
sence of inter-qubit exchange and Coulomb interactions. Each 
HQ is operated at e* = e* = ej". The states of a HQ are 
shown on the left (cf. Fig. [^; the qubit states are labeled 
by I 1 ) and j O) . The leakage states are the excited, hy¬ 
bridized state between the (2,1) and (1, 2) configurations |x) 
and the S = 3/2, Sz = 1/2 triplet state |y). The two-qubit 
states are shown on the right. The computational subspace 
I j , I , j , I j- is well separated in en¬ 

ergy from all the leakage states | q) , | /3), | 7 ), | 5), and | v) 
that are defined in Tab. Q] 


the ( 2 , 1 ) and ( 1 , 2 ) configurations at QD^ and QD 2 or 
at QD]^ and QD^. These states have lower energies than 
the state with S = 3/2 of the same (cf. Fig.j^ 1 1 ) and 
|0) have S' = 1/2, |y) has S = 3/2). As a consequence, 
all the computational states have lower energies than the 


-;rL -ttR 


~^l-j -^tv I I 

leakage states (unless <5^,(5^ < ^). Fig. 3 shows the 

energies of the single-qubit states and the le^age states 
for two identical HQs. 

A tunnel coupling between QD 2 and QDf couples the 
(2,1) and (1,2) configurations at these QDs. Similar to 
the single-qubit interactions in Sec. [IT] I phenomenologi¬ 
cally introduce a state coupling between QD 2 and QD]^: 

n'= -h\ts) (10) 

+ to It s) (S tl +h.c. 

Note that only the states with S = 1/2, Sz = 1/2 at QD 2 
and QD]^ are coupled in Eq. (10 1 . All triplets at QD]^ are 


neglected because these states have much higher energies. 

The setup can be tuned towards the (1,1,2,2) 
configuration. Only states with a singlet at 
qdR are considered, giving the four states 
{|ttSS),|ttSTo) ,|nST+),|/tST+)} with 

Sz = 1. If the ( 1 , 1 , 2 , 2 ) spin configuration is only 
virtually excited, then can be eliminated using 
Schrieffer-Wolff perturbation theory.l^sHIll Only the effec¬ 
tive interactions on the qubit subspace are summarized 
in the following; Appx. [^ gives a detailed summary of 
their derivations. 

There is an effective two-qubit interaction 


£2Q = 


7^ 

1 


( 11 ) 


'two qubits 
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Table I. 
with S 

and I vq) 

\o)AP): 


Relevant two-qubit states with total Sz = 1 for two HQs. All the states of a three-electron DQD 

= 1/2, Sz = ±1/2 have contributions in the (2,1) and (1,2) configurations. The computational subspace 

, I , I , I I is energetically separated from all the leakage states. I use the abbreviations |hi) 

for the basis states in 5 = 1/2, Sz = —1/2 that are obtained from | l) and | O) by flipping all the spins. The states 
17), I 5), and I u) label the leakage states with the energies ^ ± 5^, — ^ ± 5^, 5^ ± and 5^ ± 5*^. 


state 


qL^R\ 

0^0^' 


(S^,Sz) (S^,s?) energy 


(14) 


i/i 

V2’ 2/ 


f Aq ± Aq j /2 
/2 

^ —Aq ± Aq J /2 
— (^2\o ± Aq ^ /2 


l«l) = 1 '') [v 5 4T+) ±^|tTg)] 

(14) 

(14) 1 

(4 

102) =|4) ltT+) 


4>-|) 

(14) J 


4i) = [y4QT+) ±yi|tTg)] 

44) 

(14) 1 


1/32) =14) ltT+) 


(4-1) 

(14) J 


|7i) = [y|UT+) ±^|tTo)' 

ll») 

f3 1 \ 

\2^ 2 J 

44) 1 

1-4 ± ^ 

I 72 ) =|tT+) |4) 


(3 S'! 
V2’ 2 ; 

(4-1) J 


14) = [y||iT+) ±^|tTg)' 

|o^) 

(3 1\ 

V2 ’ 2/ 

(44 1 

[4 - ^ 

14) = |tT+) |4) 


(14) 

(- --) J 

V 2 ’ 2 / 


4i) = [^|;t+) ±yiitTg)' 

[v^UT+) ± 

^ItTo)] (1,4 

(3 I'j > 

^ 2 ’ 2 / 


14 = /||tT_) ±y|UTo) 

ltT+) 

(3 _1\ 
V2’ 2/ 

(3 S'! 

\ 2 ’ 2 / 

»4±4 

+ 

1 

<- 

II 

4To)] 

(3 S'! 

V2 ’ 2/ 

(3 _1) 

\2 ’ 2 / > 



with X£ 


Jo ± 


)(3-llc' 


108 


Ji and Zg — 


'=l(1-4) 

6 
(2 

-4=- are the effective coupling constants 

-^( 1 , 1 , 2 , 2 )—-^( 1 , 2 , 1 , 2 ) ^ 

that are determined by the energy difference £^(17,2,2) — 

£(1,27.2) between (1,1, 2, 2) and (1,2,1,2). cp [cr] is the 

abbreviation for cos ( 4 / 2 ) [cos ( 4 / 2 ) ]• Note that there 

is additionally a single-qubit Hamiltonian 


V JqJi- Jo — 




and Ji 


PlQ —ZgC 


\^L L 




( 12 ) 


leakage states, these interactions introduce minor effects 
for the qubit’s time evolution, and they can be neglected. 


IV. UNIVERSAL QUBIT CONTROL 


I discuss gate operations for the parameters (A) Aq « 

_ 

Aq and (B) Aq Aq . Note that the HQs are always 
operated at their sweet spots and e^. 


with = 

cl{3+llcl) 

108 


cl(i+4) 

6 


V JoJi, Xg 


1+ct 


Jo - 


Ji, and = 


i-ci 


Jo + 


Besides Eq. (Ill and Eq. ( 12 ), the effective Hamilto- 


4 ( 3 - 114 ) 

108 


Jl- 


nian also contains contributions for the leakage states and 
between leakage and computational states. All these con¬ 
tributions will be included in the numerical calculations 
that contain inter-qubit exchange interactions. Assuming 
large energy separations between the computational and 


A. Nearly identical qubits Ag ~ Ag 


I consider the Hamiltonian: 


-Ha = 


^0 


“0 _R I 

— CT,, + Xa^a^ . 


( 13 ) 


If HQ'^ and HQ^ are in resonance, = Aq — Aq, then 
HQ'^ and HQ^ entangle under static time evolutions. In 
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time iSWAP [ns] 



Figure 4. Gate operations for two HQs with similar eigen- 
frequencies Aq ~ that use the static Coulomb interaction 
according to Eq. ([^, with Xc/h ~ 1.3 MHz, and the tunable 
exchange interactions from Eq. 0 -(|12[). A small detuning 
Aq — Aq S> Xc is sufficient to neglect the interactions be¬ 
tween HQ'" and HQ^ such that the qubits evolve indepen¬ 
dently. ID shows the infidelity, according to Eq. (All, of a 
static time evolution with the given parameters and A = Ac, 
which is compared to the ideal time evolution with A = 0 
after the time t = hjiAXc). Adding a small tunnel cou¬ 
pling between QD 2 and QD^ and tuning the spin conhgura- 
tion towards (1,1, 2, 2) brings HQ'" into resonance with HQ^ 
(cf. Eq. (121 , ZY = A^ — Aq). The HQs entangle after 
t = h/(4X) for A = Ac -f Ac. iSWAP shows the infidelity, 
according to Eq. of the entangling operation with the 

given parameters, which is compared to an ideal entangling 
operation according to Eq. ( |13| |. The qubits can be oper¬ 
ated independently for |a^ — Aq j /h > 20 MHz. When the 

inter-qubit exchange interaction increases, HQ'" and HQ''" are 
brought into resonance to construct high-fidelity entangling 
operations. Note that also leakage errors increase with the 
inter-qubit exchange interactions, as discussed in the text. 


this case, the two-qubit interaction is 


A 


o/rwa _ 

— W Wx j 


(14) 


in the rotating frame with ^ (cr^ + cr^) for A. 

After the time t = the iSWAP gate is created, which 
is maximally entangling. l^^^^^ On the other hand, the two- 


qubit interaction in Eq. (131 can be neglected for Ag 
yf Aq A and all single-qubit gates can be realized 
with the methods that were described in Sec. El 

Altogether, universal qubit control requires rapid con¬ 
trol over Z^. For Z^ = 0, all the single-qubit gates can 
be realized with resonant drivings of the qubits at their 
eigenfrequencies. For Z^ = Ag — Ag, the HQs entangle 
under static time evolutions. 

In the configurations of Fig. A in Eq. (131 has contri¬ 
butions from the capacitive couplings between the charge 


configurations Ac (cf. Sec. Ill AI and from the ex chang e 
interaction Af between QD 2 and QD^ (cf. Sec. IIIBl: 
A = Xc + Xg. The exchange interaction causes also a 
frequency shift of HQ'" [Z^ in Eq. @1 which can be con¬ 
trolled electrically. Note that the Coulomb interaction re¬ 
mains constant during the qubit manipulations because 
the HQs are always operated at their sweet spots. The 
single-qubit interactions from Eq. (|^ are neglected be¬ 
cause they only modify the positions of the anticrossings 

and e^. 

Ag and Ag should have similar magnitudes such that 
only small inter-qubit exchange interactions are needed 
to bring HQ'" and HQ'^ into resonance. Still, Ag and Ag 
are sufficiently distinct for independent single-qubit con¬ 
trol. The capacitive coupling Xc between the HQs must 
remain much smaller than in the naive approximations 
in Sec. |HI A] Fig. shows simulations with k = 0.04 /reV 
and all the parameters defined earlier. For the Coulomb 
interactions, only the two-qubit interaction Xc is taken 
into account according to Eq. The exchange inter¬ 
action contains all the contributions from Eq. 0 and 
Eq. (@. 

If the exchange interactions between HQ'" and HQ'^ 
\Zg in Eq. (12) and Z'" in Eq. @1 cancel the detun¬ 
ing between the HQs, then the entangling gate is real¬ 
ized. Simultaneously with the single-qubit energy shift, 
also Xc and from Eq. (12) increase. These interac¬ 
tions cause systematic gate errors. Fig. shows that 
the gate errors of the iSWAP below 1% can be real¬ 
ized for Ag — Ag /h < 80 MHz. Leakage errors are 

also caused by the exchange couplings between HQ'" and 
HQ'^, and these leakage events significantly contribute 
to the gate infidelities. Fig. [^extracts the leakage errors 
during an entangling gate with finite tunnel couplings 
between HQ'" and HQ'^. The leakage probability is ex¬ 
tracted from the time evolution U by taking the norm of 
the matrix \U-p^-p^\^ between the states from the compu¬ 
tational subspace Vc and the leakage subspace Vc- 

If Z'" is reduced, then HQ'" and HQ'^ decouple. Fig. 
compares the time evolution of Eq. (131 with the time 
evolution of 


0/ideal 
ft A 


^0 L I ^0 R 

H-■ 

2 2 


The differences rise significantly if 


—r-L 

^0 ~ ^0 


but the infidelities stay below 1% for 
20 MHz. 


—r—L “a"?!' 

^0 ~ ^0 


(15) 

decreases; 
/h > 


B. Distinct qubits Aq > Aq 


_ 

Highly detuned HQs (e.g. Ag ^ Ag ) can be operated 
exclusively with microwave signals. The o-^cr^ interac¬ 
tion between HQ'" and HQ'^ can be neglected without 
drivings, and if both qubits are driven with their own 
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time CPHASE [ns] 

40 60 80 100 120 



Figure 5. Resonant manipulations of two HQs with large 
detunings between their eigenfrequencies Aq /h and Aq /h oi 
a few GHz. The qubits interact via the Coulomb interaction, 
according to Eq. with Xc/h = 0.13 GHz. A resonant 
drive of HQ'^ at its eigenfrequency realizes a state inversion 
of HQ*", while HQ^ evolves trivially. The NOT*" gate is re- 


(Aot/h) 


alized after the time t = h/(4:A) for a drive ^cos 

with A/h = 0.1 GHz. The time evolution is compered to an 
ideal NOT*^ with Xc = 0, showing that the gate infidelity, 
according to Eq. (All, is lower than 1%. A similar drive 
of HQ'" with the frequency Aq /h entangles HQ'^ and HQ^. 
The resulting gate operation is equivalent to a CPHASE gate 
up to single-qubit gates. In the simulations, the gate times 
of the entangling gates are optimized numerically to cancel 
high-frequency fluctuations. The gate errors are quantihed 
by the deviations of the Makhlin invariants from their ideal 

values (|Gi| -b IGa - ll).!^ 
high-hdelity CPHASE gates. 


Aq — A^ jh > 2 GHz permits 


resonance frequencies. In the cross-resonance protocol, 
one qub it is driven at the eigenfrequency of the other 
qubit 


■e.g., 


■Hb = -b vAcOS 




a'' 


•? -b X<7^(J^. 


(16) 


Transforming Eq. (161 to the rotating frame with 


{a\ -b cr^) gives after a rotating wave approximation 
A T T 


-^L 7 ^^ 

o/rwa ^0 ^0 _L , 

- -2-+ 2 


<7^ + 


/ L R, I L R'\ 
\^x j ■ 


(17) 


Another rotating wave approximation gives the effective 
interaction 



X [cos {d) + sin (r^) a^] 


in the rotating frame with 
abbreviations cos {A) = 


I use the 


^(a^-a^Y+a^ 


and sin {'d) = 


in Eq. (I 8 I. A drive for the time t = 

(a^o-^oY+^^ 

^ I ^ ^ ^ ^ a entangles the HQs because the time 




Y{a^-a^Y+a^ 


evolution of Eq. (181 is equivalent to a CPHASE gate.*^ 

This entangling operation only requires the in¬ 

teraction, which can be obtained from the Coulomb in¬ 
teraction between the HQs [see Eq. (|^]. The numerical 
simulations use k = 4 /ieV, which is smaller than the 
suggested value from the naive estimates in Sec. [HI A| 
Because the inter-qubit exchange interactions are not 
needed, which would require QDs in close vicinity, one 
can reduce k in Eq. ([^ by sufficiently separating the 
DQDs. In this case, there are no leakage errors (see 
Sec. HIBI because the electron transfer between the QDs 
is forbidden. I simulate the time evolutions according to 
Eq. (16) without any rotating wave approximations. Be¬ 
sides the dominant time evolution from Eq. (181, there 
are rapidly oscillating terms that were neglected in the 
rotating wave approximations in Eq. ( |17| ) and Eq. (18). 

Fig. [5] shows the errors of entangling operations for 
X = 0.55 /reV, driving amplitudes A/h — 0.1 GHz, 
and the HQ parameters defined earlier. To realize 
high-fidelity CPHASE operations, the deviations of the 
Makhlin invariants from Gi = 0 and G 2 = 1 are mini¬ 
mized. This routine cancels artifacts that are neglected 
in the rotating wave approximations of Eq. (161 and 
Eq. (HT 


If the HQs are detuned by several GHz, then 
high-fidelity entangling operations can be realized. Fig.[^ 
also shows that high-fidelity single-qubit operations can 
be realized in the same setup if the HQs’ eigenfrequencies 
are detuned by a few GHz. 


V. NOISE DISCUSSION 

This section shows that noise from fluctuating mag¬ 
netic and electric fields only causes small errors for HQs. 
In the following, a HQ is analyzed in the eigenbasis 
{1 1 ) , I 0 ) } according to Eq. (§ because it is operated 
at its sweet spot e*. 


A. Hyperflne Interactions 

Fluctuating local magnetic fields were identified as a 
natural problem for spin qubits.l^ The nuclear magnetic 
fields of the host’s nuclei couples via the contact hyperfine 
interaction to the spin of an electron.l^^lt^ For localized 
electrons, the contact hyperfine interaction can be de¬ 
scribed by slowly-fluctuating local magnetic fields at the 
QDs: and 53^°^. 

For the HQ, a noise term -b 

describes the fluctuating magnetic fields. 
























where = (cr^ 

tors for the electrons at QD^. 
gives the contribution 

jSB 


QD 
O' z 


) are the Pauli opera- 
For the HQ, this term 


'^{|T),|o)} - 


(19) 


with Sb = 


2 


_ 3+cos^(f) ^^QD, ^ 8cos^(f) ^^QD, 


and (Tj, = |l) (0| + |0) (l|. 

The nuclear magnetic field can be treated as static dur¬ 
ing a single gate operation, but it fiuctuates between suc¬ 
cessive measurements. Typical magnitudes of are 

5 mT for GaAs QDs and 100 /iT for Si QDs. The as¬ 
sociate frequencies are /h « 30 MHz for GaAs 

QDs and gg.BSB^^/h « 3 MHz for Si QDs, while Aq/Zi 
reaches several GHz. Therefore, one can treat Eq. (19) 
as a small perturbation to Eq. ([^. 


B. Charge Noise 


Fluctuating electric fi elds a re caused by the filling and 
unfilling of charge trap^^lllll or by fluctuating gate po¬ 
tentials. One can describe the dominant contribution for 
a DQD qubit by introducing an uncertainty in the de¬ 
tuning parameter which causes the noise term 




6^ 


{|i))o)} 2 ((52 +A?) 


6e 


( 20 ) 


for the HQ at the sweet spot e*, with ax = |l) (0| + 

| 0 > (T|- 

Gharge fluctuation are slow compared to the gate 
times. The rms values of 6e reach a few geY (1 geV/h « 
0.24 GHz) in GaAd^^^^^I^ and Si.^^^The factor 


!(5MAf) 

in Eq. (|20| gives an algebraic suppression to charge noise 


compared to a charge qubit. Note that this suppression 
is much smaller than for spin qubits.l^The parameters of 
the previous analysis suggest magnitudes of the fluctua¬ 
tions . / 1 /ieV//i « 76 MHz. This number it still by 


2{s^+Al) - 

two orders of magnitudes smaller than Aq/H « 4.5 GHz. 


VI. SUMMARY AND CONCLUSION 

This article has analyzed quantum computation for the 
HQ. The HQ is an exchange-only qubit,!^ where three 
electrons are confined at a DQD.l^ There are two avoided 
level crossings between states in (2,1) and (1, 2) in the 
transition region of these charge configurations. Con¬ 
trolled transfers through the avoided crossing have re¬ 
alized single-qubit gates experimentally,li2HII] while the 
possibility of resonant gates mot ivates the search for op¬ 
timal operation points.^^lEIEolMl j introduce such an opti¬ 
mal operation point at a sweet spot that is exceptionally 
noise insensitive. All the gate operations can be done 
close to this sweet spot. 


This paper has derived an effective qubit description 
at the sweet spot. Two HQs can be coupled via Coulomb 
and exchange interactions. If both HQs are operated at 
their sweet spots, then the Coulomb interaction stays 
constant but the inter-qubit exchange interaction can be 
manipulated quickly. In a first approach, two qubits of 
similar eigenfrequencies are analyzed. The exchange in¬ 
teraction can bring the qubits in and out of resonance. 
Two qubits of identical eigenfrequencies entangle under 
static time evolutions, but two distinct qubits evolve in¬ 
dependently. In a second approach, two highly distinct 
qubits are analyzed. The Coulomb interaction stays con¬ 
stant, and only resonant electric signals are needed to 
realize single-qubit and two-qubit gates. 

This paper has simulated HQs in Si and GaAs 
with realistic parameters that are extracted from 
experiments Most critically, the gate operations re¬ 
quire Coulomb interactions between HQs that are small 
compared to the naive estimates for the setup. Espe¬ 
cially, if inter-qubit exchange interactions are used for 
the entangling operations, then the inter-qubit Coulomb 
couplings should be small. This paper has discussed how 
weak inter-qubit Coulomb couplings can be realized. It 
should be possible to realize such weak inter-qubit cou¬ 
plings with a careful design of the QD layout. 

Fault-tolerant quantum computation requires high- 
fidehty quantum gates with error probabilities below 
l^ |7i|72| ]y[-y giniulations showed that quantum compu¬ 
tation with this infidelities is possible if multi-qubit ar¬ 
rangements of HQ with the described parameters can be 
fabricated. Nuclear spin noise and charge noise are less 
critical. Additionally, these statistical errors can be re¬ 
duced with refocusing proto cols w ith similar approaches 
as for DQDpHZH and TQDs.ia^lZQI A reduction of the nu¬ 
clear spin fluctu ation can be realized by preparing the 
nuclear spin batlP^^or with QD materials that contain 
nuclei of zero spin d^^ES I jf the need arises. 


A recent experiment also realized resonant single-qubit 
gates for the HQ.l^In this case, the HQ is operated deep 
in (1, 2) which reduces the influence of charge noise. Res¬ 
onant single-qubit gates are possible when the transition 
region to (2,1) is approached, while one still stays away 
from the anticrossings. Two-qubit gates with exchange 
interactions cannot be used to entangle HQs in Ref. m 
in the same way as in the study of this paper because 
many leakage states are degenerate with the computa¬ 
tional states. In Ref. izni, the S' = 3/2, S 2 = 1/2 state is 
nearly degenerate with the qubit states (cf. Fig.|^. The 
Goulomb interaction can still be used to entangle HQs, 
similar to Sec. |HI A[ while the interaction Hamiltonian is 
(X a^a^. Universal qubit control likely requires DG con¬ 
trol of the interaction Hamiltonian. Either two HQs are 
isolated from each other, or they ar e coup led via a^af- 
(in contrast to the approach in Sec. IV B| where single¬ 
qubit and two-qubit gates are possible without changing 
the operation points of the HQs). 

Overall, the HQ is an interesting candidate for fur¬ 
ther experimental and theoretical studies. It has many 
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characteristics of a charge qubit, especially with its fast 
operation times. On the other hand, it can be protected 
from charge noise similar to a spin qubit. The described 
sweet spot manipulations classify the HQ as a mixture of 
a charge qubit and a spin qubit that has advantages from 
both setups. This paper has shown that universal gate 
operations can be realized for the HQ at the threshold of 
quantum error correction, which should further motivate 
the search for optimal manipulation protocols of HQs. 
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Appendix A: Gate Fidelities 


The fidelity that is used in the paper should be defined 
in the following. To characterize gate fidelities, the time 
evolution U is compared to the ideal time evolution Ux- 
The state space is doubled to two identical Hilbert spaces 
R and S. The entanglement fidelitj®! 

F = tr{p^®l^( 8 ) [[/-^C/i]®} (Al) 


is a measure for the gate performance. Eq. (All com- 
pares the time evolution ^C/] of S with a reference 
system R that is unchanged, = \\jj) is a maxi¬ 
mally entangled state of the combined Hilbert space, e.g., 

I ^/>) = (1 0000) -f 1 0110) -f 1 1001) -f I nil)) /2. 


Appendix B: Effective Inter-Qubit Exchange 
Hamiltonian 


This section describes the derivation of the effec¬ 
tive interactions from Sec. IHIBI for weak tunnel cou¬ 
plings between HQ'^ and HQ^. The HQs are oper¬ 
ated at their sweet spots. The dominant Hamiltonian 
H is determined by Eq. 0 [with the approximations in 
Eq. ([^]. Eq. (10 1 defines the coupling between QD 2 
and QD^ (called 'H'). The system is tuned towards 
( 1 , 1 , 2 , 2 ), while the computational states still remain in 
the ground states configurations. The low-energy sub¬ 
space V = Vq + Vc contains the qubit states Vq = 

I l^T^^ , , 0^1^^ , 0^0^^ |. The leakage sub¬ 


space contains the states Vc = {| a) >1/5) , 1 7 ) , | <5) ^\v)} 
that are defined in Tab. |T 1 

Only states that have a singlet at 
QD5^ are included in (1,1,2,2): Q = 

{|ttSS),|ttSTo) ,|nST+),|itST+)}. A weak 
tuning towards (1,1,2,2) only occupies Q virtually. 
Eq. ( |Io| ) couples states in ( 1 , 2 , 1 , 2 ) and ( 1 , 1 , 2 , 2 ). 
Schrieffer-Wolff perturbation theory constructs an 
effective Hamiltonian on V, assuming that the coupling 
between V and Q is weak. Additionally, Q has higher 
energy than "P.mmil 

The effective interaction on V is defined in second- 
order Schrieffer-Wolff perturbation theory by 


n-p — rip 

+ 2 ^ 

^ keQ 


\Ei — E/fc 



Eq. (Bll uses the transition matrix element between the 


states |f) and |j) from 7^: T-Cp = {i\'H\j). and 


are the transition matrix elements between the 
states |f) from V and \k) from Q. Ei is the energy of 
the state | i). 

In the whole paper, it is assumed that the energy dif¬ 
ference between V and Q is large, especially much larger 
than the energy differences between states in V or be¬ 
tween states in Q. One can therefore use for all \ i) G V 
and \k) G Q: Ei — Ek ~ E'(i 2,1,2) ~ £'(1,1,2,2)- 
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